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ABSTRACT 

The long-duration gamma-ray burst (GRB) 030329, associated with supernova (SN) 
2003dh, occurred inside a star-forming dwarf galaxy at redshift z — 0.1685. The low 
redshift, and a rich set of archival Hubble Space Telescope (HST) images, makes this 
GRB well-suited for a detailed study of the stellar population in the immediate vicinity 
of the explosion. Since the lifetime of a star is directly tied to its mass, the age of the 
stellar population can be used to put constraints on the GRB and SN progenitor mass. 
From the HST images we extract the colours of the precise site from which the GRB 
originated, and find that the colours are significantly different from those of the overall 
host galaxy and the surrounding starburst environment. 

We have used spectral evolutionary models, including nebular emission, to care- 
fully constrain the age of the stellar population, and hence the progenitor, at the very 
explosion site. For instantaneous burst models we find that a population age of ^5 
Myr best matches the data, suggesting a very massive (Mzams > 50 M Q ) star as the 
progenitor, with an upper limit of 8 Myr (Mzams > 25 M Q ). For more extended star 
formation scenarios, the inferred progenitor age is in most cases still very young (< 7 
Myr, implying Mzams > 25 M Q ), with an upper limit of 20 Myr (M Z ams > 12 M Q ). 

These age estimates are an order of magnitude lower than the ages inferred from 
the overall host galaxy colours, indicating that progenitor mass estimates based on 
data for spatially unresolved GRB host galaxies will in general be of limited use. Our 
results are consistent with the collapsar scenario. 

Key words: Gamma rays: bursts - galaxies: starburst - galaxies: stellar content - 
stars: evolution GRB 030329, SN 2003dh 



1 INTRODUCTION 

It is now well established that long duration GRBs are con- 
nected to the deat hs of massive stars. The st andard pic- 
ture is reviewed by IWooslev and Blooml (120061). where the 
preva iling model is a collapsar ( MacFadvden fc Wooslevl 
1999), i.e. a massive star collapsing to a black hole and 
producing an energetic supernova. The first piece of evi- 
dence for the supernova scenario came with t he discovery of 
SN 1998bw in the error box of GRB 980425 (|Galama et al.l 
1998; iPatat et al.ll200ll ). and was later confi rmed by the 
association of GRB 030329 with SN2003dh |Hiorth et all 
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2003; IStanek et all 120031 : iMatheson etall boOS). It is gen- 
erally assumed that the progenitors to these GRBs are very 
massive stars. Not all types of core-collapse SNe produce 
GRBs, and it is thought that only the upper mass range 
of the IMF is contributing to the GRBs. Moreover, the 
modeling of the properties of the supernovae accompany- 
ing the long GRBs points to very massive progenitors (e.g. 
IWooslev and BloomflioOrj ) . In the case of SN 2003dh, which 
was associated with GRB 030329, estimates in the range of 
25 to 40 Mq has been obtained from modelling the sup er- 
nova observations (jMazzali et al.l [20031 : iDeng et al.l l2005). 

Arguments for a massive stellar progenitor for long- 
dur ation GRBs also com e from studies of the host galax- 
ies. iFruchter et alj (2006) confirmed that the GRBs occur 
in blue star-forming galaxies, and are concentrated to the 
most active star forming regions. For some nearby super- 
novae not connected to GRBs, a study of the actual progen- 
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itor star has been possible. Recent programmes have taken a 
systematic approach to this question and have revealed sev- 
eral supernova prog e nitor stars in p r e-explosion images (e.g . 
ISmartt et all |2004 iLi et all 120051 : iGal-Yam et aH 120071) . 
However, all the direct detections so far are for type II SNe 
with moderately massive progenitor stars. Furthermore, the 
progenitors of type Ib/c SNe st i ll remain to be direct ly de- 
tected (e.g. iMaund et all 120051 ; ICrockett etafl l2007t ) with 
the exception of SN 2006jc for which an outburst was ob- 
served c oincident with the SN just two years b efore its ex- 
plosion (|Pastorello et alj2007t [Foley et al.ll2007f ). For GRBs, 
which are typically more distant, we cannot detect single 
stars prior to explosion. An alternative way to probe the 
progenitor is to analyse the population of stars near the ex- 
plosion site. Such modeling of the entire host galaxies has 
revealed that, overall, the stellar populations of GRB hosts 
are dominated by young stars. However, to put constraints 
on the actual GRB progenitor requires a close look at the 
very site where the GRB was born and exploded. A pop- 
ulation stu dy of the environment of GRB 980425 was car- 
ried out bv lSollerman et ail ((2005), who concluded that this 
GRB orig inated in a s t ar mo re massive than 30 Mq. More 
recently, iThone et alj (|2007h found support for a massive 
star origin of the GRB 060505 which, however, did not have 
an associated luminous supernova. These results are consis- 
tent with the standard picture, although further and tighter 
constraints should be sought. 

To probe the immediate environment of the explo- 
sions requires nearby bursts and/or high-resolution imaging. 
In this paper we investigate archival Hubble Space Tele- 
scope (HST ) images obtained of t he host of GRB 030329 at 
z — 0.1685 (jGorosabel et alj|2005f ). The host appears to be a 
subluminous (L ~ 0.016L*), metal-poor, galaxy with a star 
formation rate of ~ 0.5 Mq yr _1 and a total stellar mass of 
M ~ 10 8 -10 9 Mq , dominated by a rather young stella r pop- 



ulation (|Gorosabel et all 120051 : ISollerman et al.ll2005l ). The 
host galaxy will be further discussed in Sect. 4.2. 

In this paper, we use the HST images to zoom in on the 
exact location of the GRB explosion. We show that useful 
constraints on the progenitor mass can be obtained by care- 
ful modeling of the colours as seen by HST. Such studies 
are complementary to e.g. modeling of the supernova com- 
ponent - and will help to constrain the properties of the 
exploding star. 

In Sect. [2] we describe the observational material used. 
In Sect. Owe describe how the optical data were analysed, 
i.e. how we subtracted the old population from the young 
star burst region. In Sect.[4]we describe the modeling of the 
stellar populations and discuss the constraints we derive on 
the progenitor. In Sect. [5] we discuss how age estimates can 
be improved and in Sect. [6] we summarise our conclusions. 

Throughout the paper we are assuming the following 
cosmology: Q M = 0.3, Qa = 0.7, H a = 72 km s _1 Mpc" 1 . 



2 OBSERVATIONS 

All the data analysed in this paper have been retrieved from 
the Hubble Space Telescope (HST) data arc hive. Some of the 
image s ha ve previously been dis cussed by ISollerman et al.l 
( 2005) and lFruchter et al l ^OOd ). 



2.1 The ACS/WFC data 

GRB 030329 was observed with the Wide Field Camera 
(WFC) of the Advanced Camera for Surveys (ACS) on sev- 
eral occasions between about two weeks after the explosion 
and more than two years subsequently. Since we are here 
interested in the properties of the host galaxy, early images 
are of little value except for pinpointing the exact GRB lo- 
cation. The images used in this paper are summarised in 
Tabled] 

The drizzled F606W images were aligned to a common 
origin We co-added these images sorted in four temporal 
bins: 'early', 'mid', 'late' , and 'very late' (see Table[T]). The 
'early' images are useful for identifying the precise location 
of the GRB, but the afterglow contaminates the light from 
the host galaxy. In the 'mid' images, the afterglow has faded 
enough that these data can be used to study the host galaxy 
far from the explosion site. The 'very late' images contain no 
detectable afterglow or supernova light, and from comparing 
'late' and 'very late' images, we conclude that 'late' images 
are safe to use also for the GRB environment. We made 
a total GRB-free image by masking out the GRB affected 
regions in the mid images and adding these to the late and 
very late images. In this way the final image has maximum 
depth in the outskirts of the host galaxy, but no afterglow 
or supernova contamination. The resulting image is shown 
in the upper left panel of Fig. [1] We used the same strategy 
to construct a GRB free image in F814W. For F435W we 
used the co-added images from a single 'late' dataset. 

Photometric zero points were taken from the ACS www- 
pagefQ and are: 25.779, 26.398, 25.501 for the F435W, 
F606W and F814W filters, respectively, all in the VEGA- 
MAG system, which will be adopted throughout this paper. 
In what follows we shall refer to magnitudes in this system 
as -B435, V606 and 7si4- All data that is presented has been 
corrected f or a foreground Galac tic reddening of E(B — V) = 
0.025 mag (|Schlegel et al.ll 19981 ). Taking the different band- 
pass of the used filters with respect to the 'standard' system 
into account, the applied corrections are (in magnitudes): 
A(B i35 ) = 0.109, A(Veo6) = 0.075, A(/ 814 ) = 0.049. 



2.2 The ACS/HRC/F250W data 

The archive also contains several observations in the F250W 
filter with the UV-sensitive High Resolution Camera (HRC) 
of ACS, which were also retreived. However, we could not 
properly align the HRC/F250W observation with those ob- 
tained with the WFC. The different apertures of the two 
cameras results in slight inconsistencies in the astrometric 
keywords. As a result, we could not be sure that the Mul- 
tidrizzle software had correctly registered the images at 
the requested pixel-level accuracy. Remaining misalignments 
could not be fixed by hand since no for- or background point- 
sources were available for the image registration. Hence, we 
did not use the F250W data in our final analysis. 



http:/ /www. stsci.edu/hst/acs 



Constraining the mass of the GRB 030329 progenitor 3 



Table 1. ACS/WFC images used (exposure time in seconds). 



Filter Date Dataset name Exp. time Class 



F435W May 24, 2004 J8IY9V010 1920 late 



F814W Nov. 12, 2003 J8IY2T050 2040 mid 

F814W May 24, 2004 J8IY9V020 2040 late 



annuli. The structure of the outskirts of the host galaxy 
is well fitted by an exponential disk with scale length of 
f/.'3 (0.9 kpc) and an inclination angle of i = 55° (assum- 
ing cos(i) = b/a, i.e. an infinitely thin disk). We created a 
synthetic galaxy image with these properties, scaled to the 
surface brightness level of the outer parts of the host galaxy 
and subtracted it from the host galaxy image. In this way 
we obtained an image of the central starburst with no con- 
tamination from the underlying galaxy. The same method 
was used to obtain images of the central starburst in the 
F435W and F814W filters. Being cleaned from the underly- 
ing disk contribution, these images are suitable for a pho- 
tometric modeling of the properties of the GRB explosion 
environment. In Fig. 1 we show the luminosity profile of the 
host and the resulting image after subtracting the underly- 
ing host galaxy. 

We then proceeded to extract colours for the slightly off- 
centre GRB explosion site. In order to check how sensitive 
the results are to the host galaxy subtraction, we increased 
the luminosity of the disk model to the upper limit ('max 
disk') allowed by the uncertainties, and also derived the 
colours without subtracting an underlying disk ('no disk'). 

This photometric method is appropriate if the stellar 
population from which the GRB originated dominates the 
light within the aperture. An alternative method is to also 
subtract the median flux in a surrounding annulus, as is 
done in conventional aperture photometry, and this method 
was used as a consistency check. 

The GRB environment photometry is presented in Ta- 
ble 2, and was extracted for an aperture radius of r = Of.'l. 
We also experimented with a smaller aperture which gave 
consistent results. For the r = O'.'l photometry, aperture cor- 
rections of 0.510, 0.461 and 0.538 mag were subtracted from 
£>435, ^606 and 7814, respectively l|Sirianni et alj|2005h . 

The scaling of the underlying disk has a minor impact 
on the colours of the GRB explosion site. The statistical un- 
certainties are estimated to be ob = 0.03, oy = 0.015, aj = 
0.035 magnitudes. The differences in Table 2 give a good idea 
about the systematic uncertainties, in which the subtraction 
(or not) of the local background has the largest impact. By 
modeling the stellar population (see Sect. 14.3)) using the full 
range of values in Table 2, we also include the systematic 
uncertainties in our analysis. 

In addition, we present in Table 3, the photometry for 
the whole galaxy. 

3.1 NICMOS data 

We were not able to clearly detect and fit the underlying 
host galaxy in the NICMOS images, but have extracted the 
flux in the GRB-region in apertures matching those of the 
WFC data. In the F160W image many pixels near the GRB 
explosion site are negative and we regard these data as un- 
reliable and will not use them further. In the F110W image 
we could extract a magnitude for the same region as used 
for the WFC data but without subtraction of the underly- 
ing galaxy (corresponding to the 'no disk' option). In this 
way we obtained Jno = 24.9 ± 0.4 mag., where we have 
adopted the definition of total NICMOS magnitudes and an 
aperture correction derived from the encircled energy distri- 
bution given in th e HST/NICMOS Instrument Handbook 
l|Barker et alj|2006l ). 



2.3 NICMOS data 

GRB 030329 was also observed with the NIC2 camera of 
the NICMOS instrument onboard HST. These data suffered 
from well known artefacts in the form of non-linear temporal 
bias drifts during the exposure^. The data were re-reduced 
using the CALNICA pipeline and the artifacts removed us- 
ing the BIASEQ and PEDSKY tasks of STSDAS. Unfor- 
tunately, the latest epochs which were the only ones free 
of GRB contamination suffered from cosmic ray persistence 
limiting the depth of these observations. Aligning the NIC2 
images with ACS/WFC was complicated due to the very 
few reference sources visible within the NIC2 field of view. 
We therefore rotated the images to the orientation of the 
WFC data and used the GRB itself, present in the early 
NIC2 images, for centring. The images potentially relevant 
for our analysis are made up of two exposures in each of 
the F110W and F160W filters, taken on May 23, 2004. The 
total exposure time in each filter is 1216 s. However, only 
the F110W image turns out to be useful for our purposes, as 
further discussed in Sect. El The adopted NIC2+F110W ze- 
ropoint is 22.95 (Vegamag) and magnitudes will hereafter 
be denoted by Jno- 

2.4 HST spectroscopic data 

GRB 030329 was observed with both STIS and ACS in spec- 
troscopic mode. It could potentially be very useful to con- 
strain the internal extinction from the Ha/H/3 ratio at the 
actual location of the burst. However, these spectr43do not 
have enough signal, or are too heavily contaminated by the 
GRB afterglow at early epochs, to allow such an analysis. 



3 PHOTOMETRIC ANALYSIS 

We examined the shape of the host galaxy in the final very 
deep F606W image. We fitted ellipses to the outer isophotes 
and derived a luminosity profile by integrating in elliptic 

2 see http:/ /ww w.stsci.edu/hst/nicmos| 

3 http://www-int.stsci.edu/~fruchter/GRB/030329/index.html 
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Figure 1. Upper left: The host galaxy of GRB 030329 in the F606W filter (HST/ACS). The field shown is 4'.' 2 X 4'.' 2 which corresponds 
to 12 X 12 kpc at the distance of GRB 030329. This image is made from data obtained at different epochs and is not contaminated by 
light from the GRB afterglow. The location of the GRB explosion is indicated by the concentric white circles (where the innermost 
circle corresponds to the used photometric aperture). Upper right: Synthetic exponential disk model fitted to the outer isophotes of 
the deep F606W image. The inclination is 55° and the scale length 0'.'31 (0.87 kpc). Lower left: Here we show the result of subtracting 
the disk model from the F606W image in the upper left panel. This has isolated a central starburst region, free of contamination from 
the underlying population described by the best fit disk model. The GRB occurred near the edge of the central starburst region, ~ 0'.'2 
from its centre. Lower right: The luminosity profile resulting from integrating the unsubtracted F606W image in the upper left panel 
in concentric ellipses (with centroid from the outer isophotes) is shown as a solid line while the fitted disk model is shown as a dashed 
line. 
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Table 2. GRB environment photometry. The upper entries are for photometry extracted in a ff.'l radius aperture without subtracting 
the local background. For the lower entries, local background has been subtracted using the median value in one pixel wide annulus 
immediately outside the photometric aperture. All magnitudes are in the HST VEGAMAG system. 
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Max disk 
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24.843 
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Best disk 


no 
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24.766 
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0.561 


No disk 


no 
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0.596 


Max disk 
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25.126 


25.156 


24.696 


-0.044 


0.480 


Best disk 
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25.114 


25.159 


24.678 


-0.029 


0.460 


No disk 


yes 


25.117 


25.116 


24.610 


0.002 


0.506 



Table 3. Integrated host galaxy photometry within the Holmberg 
radius (1.4", equivalent to 3.9 kpc) 



-B435 V&06 hl4 -B435 — Ve06 ^606 — h'S5 

23.2 22.8 22.3 0.36 0.52 



4 STELLAR POPULATION ANALYSIS 

We will now describe how we use spectral evolutionary syn- 
thesis models together with the derived photometry to esti- 
mate the age of the stellar population. Once an age estimate 
has been found, this can be converted into a lower limit on 
the zero age main sequence proge nitor mass (M z ams) using 
stellar evoluti o nary models from iMevnet et al.l (|l994l ) and 
iFagotto et alj (|l994T ). The lifetimes for these two indepen- 
dent models agree well, typically to within 5%. More recent 
calculations (jHirschi et al.ll2004l ) that include stellar rota- 
tion indicate that, at least for solar metallicity, the effects on 
the life times are insignifica nt. For hypothetical , extr emely 
rapidly rotating stars (e.g. IWooslev and Hegerl I200r3 ). the 
mass lifetime-relation could possibly be affected, but no 
quantitative estimates for such stars are yet availa ble. In 
what follows we will adopt the IMevnet et al.l ||l994) mod- 
els. The models assume single star evolution, but even if the 
GRB progenitor were part of a binary system, this would 
have no major impact on the lifetime which is dominated 
by the main sequence phase. 

4.1 Spectral synthesis 

Integrated photometry is commonly analysed using spectral 
evolutionary synthesis models. In this procedure, it is impor- 
tant to adopt a model as closely adapted to the known phys- 
ical properties of the target system as possible. As the GRB 
03032 9 host galaxy is a blu e object with strong emission- 
lines (jGorosabel et al.l I2005T) . nebular emission is likely to 
give an important contribution to the optical spectrum. 
Many publicly available models do not include nebular emis- 
sion and are therefore not suitable for analysing systems of 
this kind. The previous GRB 030329 hos t galaxy analysis by 
iGorosabel et al.l (|2005h is based on the iBruzual fc Charlotl 

|2003t ) model, which suffers from t his shortcoming. 

In the following, we adopt the Za ckrisson et al.l (J2001) 
spectral evolutionary model, which uses the phot oionisa- 
tion code Cloudy version 90.05 l|Ferland et al.lll99Sl ) to pre- 



dict the nebular continuum and emission lines. For each 
time step, the spectral energy distribution of the stel- 
lar population is used as input to Cloudy, which makes 
the computation time-consuming, but provides a realis- 
tic evolution of the nebular component. The stellar com- 
ponent, which includes pre-main sequence evolution and 
a stochastic treatment of horizontal branch morphologies 
at low me tallicities, is based o n sy nthetic stellar atmo- 
sphere s by iLeieune et alj (|l99Sf ) and IClegg fc Middlemassl 
(1987), together with stellar evolutionary tracks mainly 
from the Geneva g roup. We have previously demon strated 
|Ostlin et al.l 120031 ) that the IZackrisson et all l|200l[ ) model 
performs much better than the publicly available vers ion 
of the PEGASE.2 code (|Fioc fc Rocca-Volmerangell 19991 ) in 
reproducing the observed colours of young super star clus- 
ters. (Although PEGASE.2 has previou sly been used in the 
analysis of long-durat ion GRB host, e.g. lSokolov et ai"1l200ll ; 
ISollerman et~al]|2005l .) This is thanks to its more sophisti- 
cated treatment of the nebular component which is known to 
contribute substantially. To allow an analysis of the broad- 
band magnitudes relevant for the current investigation, the 
spectra predicted by the models are redshifted to z — 0.1685 
and convolved with the throughput curves of the filters used. 

In the lZackrisson et al.l ((2001) model, the properties of 
the stellar component are regulated by the metallicity, ini- 
tial mass function (IMF) and star formation history. We 
here assume the metallicity to be the same for the stars and 
the ionized gas, i.e. Z — Z B t avs = Z gas and have explored 
Z = 0.001,0.004 and 0.008. For IMF we have assumed a 
power- law (dN/dM tx M~ a ) throughout the 0.08-120 M 
stellar mass range, and we note that the actual slope at 
low masses (M < IMq) will not affect the results for the 
young ages here considered. Our standard assumption is a 
Salpeter (a = 2.35) IMF but we have also explored a = 1.85 
and 2.85. For simplicity, we assume the star formation rate 
(SFR) to be either instantaneous or exponentially decreasing 
over time (SFR(i) oc exp — t/r), where the e-folding decay 
rate r describes the duration of the star formation episode. 
An instantaneous burst corresponds to r ~ 0. 

The properties of the nebular component are deter- 
mined by the adopted hydrogen number density n(H), which 
is assumed to be constant within the nebula, the gas filling 
factor /, and the gas mass M gaa available for star forma- 
tion throughout the star formation episode. M gaB is used 
together with the adopted SFR(£) to determine the stellar 
Lyman continuum flux density for each time step. This be- 
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comes an important parameter, because the homology rela- 
tion betwee n plane-parallel nebulae with the same ionisation 
parameter [|Davidsor]|l977h breaks down for the spherical 

nebulae assumed here. 

The IZackrisson et al.l (2005) model assumes that the 
nebula responds instantaneously to changes in the ionizing 
flux from the stellar population. This approximation holds 
as long as the many processes in the ionized cloud occur 
on timescales shorter than the age resolution that we are 
interested in. To ensure that this is not an issue, we used 
Cloudy to verify that the timescales of the longest radia- 
tive processes in the ionised gas are shorter than the time 
resolution of the stellar population model. For the range of 
gas parameters used in this paper, all processes occur on 
timescales much lower than 1 Myr, which is quite sufficient 
for our purposes. 



4.2 The GRB 030329 host galaxy 

The host galaxy of GRB 030329 has bee n investigated in 
severa l studies using ground based data. iGorosabel et al.l 
(2005) used spectroscopy of t he host to estimate a metal- 
licity of Z = 0.004. However, (|Sollerman et al.ll2005l ) found 
that the available spectroscopy in the literature does not 
strongly constrain the metallicity of the host galaxy, and 
neither Z = 001 n or Z = 0.008 could be ruled out. 
IGorosabel et~ail (|2005h could not tightly constrain the ex- 
tinction but found that the host galaxy spectral energy 
distribution and Balm er line ratios favoured Ay ^ 0.6. 
iMatheson et al.l (|2003l 1 used the colour of the GRB after- 
glow to infer Av ~ 0.12±0.22. Extinction may vary spatially 
over a galaxy, and the GRB environment may well have a 
substantially lower extinction than that of the luminosity- 
weighted average of the host galaxy. This scenario also turns 
out to be consistent with our modelling results presented in 
Sect. 4.3. 

In Fig. [5] we compare the colours of the overall host 
galaxy (filled triangle) to those of the GRB explosion site 
(filled square and circle, indicating colours with and with- 
out subtraction of the local background, respectively). While 
similar in V(,oa — Isu, the -B435 — V606 colour measured at 
the explosion site is bluer than for the host galaxy by ~ 0.3 
magnitudes. Th e dereddening vectors relevant for the over- 
all host galaxy (Gorosabel et al .112005! ) and the GRB explo- 
sion sites (Mat heson et al.ll2003l ) are plotted as arrows. For 
the host galaxy, the two arrows included correspond to the 
extinct ion curves for the Milky Way and Small Magellanic 
Cloud (Pei 1992j), respectively. For the GRB site, only one 
arrow has been plotted, because the two corrections overlap 
at this resolution. After correction, the GRB site is bluer 
than the host galaxy in -B435 — V606 but redder in V6O6 — Isu- 
As will be demonstrated, this is consistent with the notion 
that the stellar population hosting the GRB progenitor is 
substantially younger than that of the overall host galaxy. 

The solid line in the left panel of Fig. [2] indicates 
the evolution predicted for a Z — 0.004, r = 10 8 yr, 
M gas = 10 s Mq stellar population with gas parameters 
/ = 0.001 and n(H) = 10 2 cm" 3 . After reddening correc- 
tions, the host galaxy -B435, Veoe and Is\4 data are con- 
sistent with a luminosity-weighted age of as 50 Myr. This 
corresponds to the lifetime of a ~ 8 Mq star, which is 
close to the minimum progenitor mass for a core-collapse 



SN l|Heger et al 1120031 ; lEldridee fc Toutll2004r ) . The collapsar 
scenario for long-duration GRBs requires a higher progenitor 
mass than this; a minimum mass of ~25-30 Mq is required 
in order for stars to have cores that collapse to black holes 
and to have suffic ient angular momentum to form a disk 
l|Heger et al.ll2003l ). 

Age estimates for composite stellar populations based 
on optical broadband photometry alone are typically not 
very precise. This is due to the dege neracy between 
star f ormation history and age (e.g. iGil de Paz fc Madord 
120021 ; IZackrisson et all 120051 ). and, for systems with ongo- 
ing star formation, because of uncertainties in the phys- 
ical conditions of th e ionised interstellar medium (e.g. 
IZackrisson et al]|200l[ ). These effects are demonstrated in 
the right-hand panel of Fig. where the -B435, V6O6 and Igi4 
data are compared to the predicted evolution for a model 
identical to the one in the left panel, except for a slightly 
shorter star formation episode (r = 5 x 10 7 yr) and a lower 
hydrogen density (ra(H) = 10 cm -3 ). In this case, the best- 
fitting age becomes as small as 1-2 Myr. The large difference 
between the evolution in the two panels of Fig. 2 at low ages 
is entirely due to changes in the nebular spectrum. Since we 
have little information on the exact star formation history 
of the host galaxy and its detailed gas properties, the optical 
broadband data simply do not contain sufficient information 
to constrain the age with very high precision. 

An age higher than 50 Myr is nonetheless favo u red b y 
the ground based near- IR data of IGorosabel et al.l (|2005^) . 
This is illustrated in Fig. [3] where the temporal evolution 
of various model sequences are compared to their V — H 
measurement for the GRB 030329 host galaxy. In all cases, 
the best-fitting ages lie around m 100 M yr, which is broadly 
consistent with IGorosabel et al.l ((2005). The ages inferred 
from V — H are relatively robust to changes in the assumed 
nebular gas parameters. 

While a best-fitting age of ~ 100 Myr would correspond 
to a GRB progenitor mass much lower than that favoured by 
the collapsar scenario, stars with much younger ages - and 
hence higher masses - would still be present in a population 
with on going star formation (i.e. r > 0). Hence, while the 
integrated host galaxy properties are broadly consistent with 
a collapsar, they do not provide useful constraints on the 
mass of the actual GRB progenitor. 

4.3 The GRB explosion site and the age of the 
GRB progenitor 

Starburst regions often contain young star clusters, some- 
times referred to as 'super star clusters' (SSCs) with 
absolute ma gnitudes M v < -10 jArp fc Sandagei Il985l ; 
iMeurer et aTTl 19951 ). At closer look, luminous starburst re- 
gions are often com posed of numerous SSCs of varying age 
(|Ostlin et al.ll2003l). Each such SSC m ust have formed on 
a short timescale (|McKee fc T an 2002) and they are com- 
monly assumed to be truly single stellar populations, even if 
the star formation history of a whole starburst region may 
be more complex. A nearby example is R136, the central 
cluste r of the 30 Doradus complex which iMassev fc Hunted 
( 1998) found to be well explained by a star formation episode 
of duration up to 4 Myrs and a standard Salpeter IMF. 

The restframe absolute magnitude of the GRB explo- 
sion site is Mv 606 ~ —14.8 which indicates a stellar popu- 
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Figure 2. The Veoe — r 8i4 vs. B435 — V606 colours of the integrated GRB 030329 host galaxy (filled triangle) and the GRB explosion site 
(filled square and circle, indicating the colours with and without subtr action of t he local background). The arrows i ndicate the dereddening 
vectors for Milky Way and Small Magellanic Cloud extinction curves ( IFeill 1992T) . given the extinction estimates of lGorosabel et al,l (2005) 
for the overall galaxy and Mathcso n et al.l ||2003|) for the explosion site. At this resolution, the two arrows overlap in the latter case. 
The smaller cross in the upper left corner indicates the 1 a statistical uncertainties of the host galaxy colours, whereas the larger cross 
indicates the uncertainties for the GRB explosion site. Left: The solid line represents the spectral evolution predicted for a burst-like 
stellar population at 2 = 0.1685 with Z = 0.004, r = 10 8 yr, Af gas = 10* M Q , / = 0.001 and n(H) = 10 2 cm -3 . The dots along this line 
indicate ages of 1 Myr, 10 Myr, 100 Myr and 14 Gyr. After reddening correction, the host galaxy is consistent with an age of ~ 50 Myr. 
Right: The dashed line represents the predicted evolution for a population with r = 5 X 10 7 yr and n(H) = 10 cm -3 , but parameter 
values otherwise identical to those in the left hand panel. In this case, the best-fitting age becomes only 1-2 Myr, demonstrating the 
degeneracies involved in determining accurate ages for young stellar populations from optical broadband photometry only. 




Figure 4. The V606 — r Sl4 vs. B435 — V6O6 colours of the GRB 030329 explosion environment (as in Fig. 2, filled square and circle represent 
colours derived with and without subtraction of the local background) compared to various model populations (lines). The extinction 
corrections are the same as in Fig. [2] Left: Stellar populations with exponentially decreasing star formation rates of r = 5 X 10 6 (thin 
dashed), 10 7 (thin solid) and 2 X 10 7 yr (thick dashed). All scenarios assume Z = 0.004, M gas = 10 7 M Q , / = 0.001 and 71(H) = 10 2 
cm -3 , leading to a best-fitting age of less than 10 Myr. Right: A stellar population with Z = 0.001, but with parameters otherwise 
identical to the r = 10 7 yr scenario from the left hand panel. In this case, the best-fitting age becomes 15 Myr. 



lation mass of ~ 10 6 Mq (the exact value depends on the 
actual age and reddening). Our extraction radius (0.1") cor- 
responds to 0.28 kpc and the GRB exploded at the eastern 
edge of the highest surface brightness region in the starburst 
(see Fig. 1 and also ?). Hence it may have occurred in a sin- 
gle very luminous SSC, or in a complex of several SSCs of 
more 'normal' luminosity. 

If the GRB occurred in a single SSC, then the star 



formation time scale (r) should be close to instantaneous 
(or have a maximum r of a few Myr), after which negative 
feedback associated with photo-ionisation and stellar winds 
would have quenched further star formation. However, our 
aperture may include more than one SSC in projection, and 
their different ages could then mimic a temporally extended 
star formation event. In the nearby blue compact galaxy 
ESO 338-04, the central 2 kpc starburst is composed of more 
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Figure 3. The temporal evolution in V — H pre dicted by various 
model s compared to the ground-based data by Gorosabcl ct al. 
(2005) for the host galaxy (thick solid line, with grey region in- 
dicating 1(7 error bars). The observed colour has been corrected 
for Galactic and intrin sic extinction - the latter assuming the 
Gorosabcl ct al ] l|2005l l extinction estimates and a Small Mag- 
ellanic Cloud extinction curve. The different model predictions 
correspond to star formation histories r = 5 X 10 7 , 10 s , 2 X 10 8 , 
3 X 10 s , 10 9 , 10 10 yr (from top to bottom at an age of ~ 5 X 10 8 
yr). A population with Z = 0.004, M gas = 10 s A/ Q , / = 0.01 and 
ri(H) = 10 2 , cm -3 has been assumed throughout. Regardless of 
model scenario, the best-fitting age lies around ~ 100 Myr. 



than 50 young SSCs with ages ranging from new-born to 
about 40 Myr (|6stlin et alj|2003l ). If the situation is similar 
in the host galaxy of GRB 030329, then the effective star- 
formation time-scale (r) in our aperture could be extended 
up to 10-20 Myr. This uncertainty on r is the most serious 
limitation for our method and is directly tied to the spatial 
resolution and distance of the target. Below we shall discuss 
the constraints on the age of the GRB progenitor for various 
values of r and other parameters. 

The left panel of Fig. U shows that the GRB explosion 
site can be well-fitted by a Z = 0.004, M gas = 10 7 M© pop- 
ulation with a star formation time scale around r = 10 7 yr. 
Contrary to the case for the overall host galaxy, the age of 
the GRB environment appears to be very young - well be- 
low 10 Myr (corresponding to Mzams > 21 Mq). Just as in 
the case of the host galaxy, the exact age is uncertain be- 
cause of the aforementioned effects related to star formation 
history and properties of the ionized interstellar medium. In 
Fig. 21 three different star formation histories are plotted 
(r = 5 x 10 6 , 10 7 and 2 x 10 7 yr). All of these scenarios 
provide reasonable fits to the observed colours of the GRB 
site, but at slightly different ages (2-5 Myr). 

The variations in the colours predicted by these scenar- 
ios at such early stages of evolution come from the different 
Lyman continuum flux densities that these stellar popula- 
tions produce. While the best-fitting ages for the GRB en- 
vironment are systematically lower than those for the over- 
all host galaxy, it is still possible to find a small number of 
model parameter combinations that would allow higher than 
10 Myr. This is demonstrated in the right panel of Fig. [4] 
where the metallicity of the r = 10 7 yr scenario from the 



left panel has been changed to Z — 0.001, giving an age of 
14 Myr. This is actually the highest best-fitting age found 
when fitting the -B435, V6O6, I»i4 and Jno magnitudes of the 
GRB environment to the 100 sets of model parameter val- 
ues explored. It corresponds to a GRB progenitor mass of 
Mzams ^ 15 M© . We once again stress, that because of the 
adopted r > star formation history, the GRB progenitor 
may in this case have a lifetime anywhere in the range from 
3 to 14 Myr, i. e . fully consistent with the collapsar scenario 
of iHeeer et all (l2003h . 

As shown in Sect. 4.2, near-IR data can be very useful 
for improving age estimates based solely on optical photom- 
etry. In Fig. [5] we compare the V^oa — Jno colours of the 
GRB environment with the evolution predicted for differ- 
ent values of r. Colours with (thick solid line) and without 
(thick dashed line) subtraction of the local background are 
shown. Due to the large photometric uncertainties, a best- 
fitting age of anywhere from to 20 Myr can be derived 
from this colour. The quality of the Jno data is obviously 
too poor to substantially improve the age estimates. 

In all, 100 models were fitted to the data in Table 2. 
The variations in colours between the different photometry 
methods give a good representation of the systematic un- 
certainties, whereas the statistical uncertainties are always 
small, with the exception of the Jno-band. We compare the 
various photometries (thus spanning the range of systematic 
uncertainties) to the models and select those solutions that 
match the data to within 3 ct of the statistical uncertainties. 
The extinction has been a free parameter in the fit, but all 
fits within 3 ct of the photometric uncertainties consistently 
produced small values: E{B^ — Veoe) < 0.1 mag. This is 
consistent with the va lue derived from the optical afterglow 
l|Matheson et alj|2003h . which is reassuring. 

For the 40 instantaneous burst models explored, we find 
that all models that provide good fits (within 3 ct of the 
photometric uncertainties) indicate a very young progenitor: 
the best fit is always less than 5 Myr (implying Mzams S3 47 
Mq) and the 3ct upper limit is 8 Myr (Af Z AMS > 26 M Q ). 

Following the discussion above, we now also consider 
60 additional models with more extended star formation 
episodes, r ^ 50 Myr. The best fitting ages are in gen- 
eral still very young with the vast majority at ^ 5 Myr. 
However a small number (1 to 3, depending on which pho- 
tometry is used) of models give decent fits for ages 7-14 Myr 
(Mzams > 30 to 15 Mq), with 3ct upper limits of 20 Myr 
(Mzams > 12 M©). One of these models has already been 
presented in Fig. 4 (right panel). 

Of all the 100 models with r < 50 Myr, 30-50% (de- 
pending on which photometry is used) provide best fit so- 
lutions that are within the 3 ct statistical limits. On average 
the photometry that provides the best fit to the models is 
the 'Best disk' without subtraction of the local background. 
The vast majority of models that produce good fits give 
very young ages (^5 Myr). However, since we cannot rule 
out the physical parameters for those few that give higher 
ages, we are left with a rather loose mass constraint for the 
progenitor. 

We note that if we would have adopted the rather com- 
mon approach of simply assuming an instantaneous burst 
with fixed metallicity and standard IMF and gas param- 
eters, our constraints would have appeared as very tight. 
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Figure 5. The temporal evolution in V606 — Jno predicted by 
various model scenarios compared to the colour of the GRB explo- 
sion site with (thick solid line) and without (thin dashed line) sub- 
traction of the local background. The open symbol in the upper 
left indicates the statistical 1 a uncertainty. The observed colours 
have been correcte d for Galactic and intr insic extinction - the 
latter assuming the lMatheson et al. I d2003h estimate and a Small 
Magellanic Cloud extinction curve. The difference between Milky 
Way, Small and Large Magellanic Cloud extinction curves (IPeil 
1 19921) only amount to an uncertainty of 0.02 magnitudes in this 
colour. The different model predictions correspond to star forma- 
tion histories r = 5 X 10 6 (solid line), 10 7 (dashed) and 2 X 10 7 yr 
(dash-dotted). A population with Z = 0.004, A/ gas = 10 7 M Q , 
/ = 0.001 and n(H) = 10 2 , cm -3 has been assumed through- 
out. Due to the uncertainty in the subtraction of the local back- 
ground, a best-fitting age of anywhere from to 20 Myr can be 
derived, and within the error bars, ages of up to 100 Myr are in 
fact allowed. Such high ages are however inconsistent with the 
#435 - Veo6 and V606 - r 8i4 colours. 

Instead, we have conservatively explored a large parameter 
space in order to really test the power of this method. 



5 OUTLOOK 

As we have shown, the GRB 030329 environment is best fit- 
ted by an age ^ 5 Myr, which is an order of magnitude 
lower than the age inferred for the overall host galaxy. This 
implies that progenitor mass estimates based on spectral or 
photometric data for poorly resolved host galaxies may be of 
limited value. Even with the high-resolution data analysed 
here, the uncertainties in the exact star formation history 
and gas properties of the GRB environment impose severe 
limits on the precision with which the mass of the GRB 
progenitor can be derived from broadband photometry. 

An additional source of uncertainty stems from the in- 
trinsic accuracy of the spectral evolutionary model used. 
Currently, the only way to assess such uncertainties would be 
to analyze the data using several different models in parallel. 
As we have shown that the treatment of nebular emission 
is very important for the age derived for the GRB 030329 
envi ronment, only models including a nebular component 
(e.g. iMov et al.ll200ll ; lAnders fc Fritze-v. AlvensleberJliooll ; 



iMagris et ail 120031 ; iDopita et all l2005t ) should be used in 
such a comparison. 

To better constrain the GRB progenitor masses more 
diagnostics are required. As already shown in Fig. 5, bet- 
ter quality near-IR data would be quite useful. This could 
be obtained with NICMOS, the next generation HST opti- 
cal/IR imager WFC3, or with ground based adaptive optics 
systems. It would also be very useful to complement the op- 
tical HST data with a V-like filter that does not transmit 
the [Oiii]aa4959,5007 emission line, since this would make the 
analysis less sensitive to assumed nebular gas parameters. 

The equivalent width of the Ha emission line, EW(Ha), 
measured through spectroscopy or narrow-band photometry, 
may also provide useful constraints. In Fig. [6] we pres e nt the 
EW(Ha) evolution predicted by the lZackrisson et alj (|200ll ) 
model during the first 50 Myr of a burst-like stellar popula- 
tion for three different star formation histories. For a given 
EW(Hcv) measurement, the youngest age is inferred by the 
shortest star formation episode (here represented by an in- 
stantaneous burst), whereas the oldest age is inferred from 
the longest period of star formation considered (r — 2 x 10 7 
yr). Even without detailed knowledge of the star forma- 
tion history, interesting age constraints can be derived from 
EW(Hq). If, for instance, EW(Ha)~ 1000 A, the allowed 
age range would be 3-8 Myr, corresponding to a progenitor 
star in mass the range 25-120 M Q . For EW(Ha)Ri 100 A, a 
lower limit of 8 Myr could be imposed, indicating a progen- 
itor mass below 25 Mq. Moreover, the EW(Hq) of a young 
stellar population is almost completely insens i tive t o many 
of the parameters used in the lZackrisson et al.1 (|200ll ) model 
to regulate the properties of the ionized gas. 

While being a powerful age indicator, the Ha equivalent 
width comes with a number of potential caveats that require 
consideration: While insensitive to the gaseous metallicity, 
there is a certain dependence on the stellar metallicity. This 
is demonstrated in Fig. [6] where we contrast the EW(Ha) 
evolution at Z stals = Z gas = 0.001, 0.004 and 0.008. With- 
out more detailed metallicity information than this, the age 
constraints derived from the EW(Ha)« 1000 A and 100 A 
measurements would be weakened to 3-20 Myr and ^ 7 
Myr, respectively. Emission-line diagnostics formally probe 
the gaseous metallicity, which may differ somewhat from the 
stellar counterpart, and this will be a lingering source of er- 
ror in the conversion from EW(Hq) to population age. 

Selective extinction, i.e. the possibility that the emis- 
sion line and the continuum originate from different spatia l 
regions with different dust content (e.g. lCalzetti et al"1ll994h . 
represents another potential complication. 

This investigation provides support for the origin of long 
durat ion GRBs in very massive stars. In ISollerman et al.l 
(2005) wc found a similar conclusion for SN 1998bw 
(GRB 980425). The parent population of core collapse su- 
pernovae is expected to be biased towards ages correspond- 
ing to the lifetime of stars with mass close to the lower limit 
for SNe (~ 8Mq), since for a normal IMF these are more 
common than very massive stars. Already the finding that 
two of the still few known SNe coincident with long duration 
GRBs seem to have higher masses suggests that the parent 
populations may be different. This is consistent with studies 
of t he spatial distributio n of GRBs and SNe within galax- 
ies ijFruchter et al.ll2006l ). It also agrees with the modeling 
of the SNe associated with GRBs, which typically indicate 
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Figure 6. The temporal evolution of the Ha emission-line equiv- 
alent width predicted for stellar populations with burst-like star 
formation histories. The thick solid lines correspond to Z = 0.004 
populations with exponentially decaying star formation rates of 
t S3 (i.e. instantaneous burst), 5 X 10 6 and 2 X 10 7 yr (from bot- 
tom to top at an age of 10 7 yr). The two thin dashed (dash-dotted) 
lines represent the corresponding equivalent width predictions for 
t ks and t = 2 X 10 7 yr when instead Z = 0.001 (Z = 0.008). 



very massive progenitors (e.g. M azzali et al.ll2007l Table 1). 
Admittedly more examples would be needed since for indi- 
vidual cases there is always a risk that photometry is af- 
fected by limited spatial resolution and chance projections 
such that the progenitor had a different age than its sur- 
roundings. By imaging the most nearby GRB host galaxies 
with HST-like resolution, population studies of the explo- 
sion environment can provide statistical constraints on GRB 
progenitor masses. The method could also be used for su- 
pernovae that are too distant for single stars in their host 
galaxies to be resolved (D > 20 Mpc) and has been used 
for constraining the progenitors of the type I IP SN 2004dj 
|Mafz-Apellaniz et al]|2004l; I Wang et al.lfSooBI ') and the type 
Ic SN 2007gr (|Crockett et al.l l2008) 



6 SUMMARY 

We have used archival HST images to derive photometry 
for the explosion site of GRB 030329. This has been used 
together with spectral evolutionary synthesis models, in- 
cluding nebular emission and a wide variety of parameter 
settings, to constrain the age of the stellar population. Tak- 
ing this age to equal the lifetime of the GRB progenitor star, 
we can use the tight relationship between initial mass and 
lifetime to weigh the progenitor. 

Studying the photometry at the resolution allowed by 
HST, the explosion site has an inferred age of ^ 5 Myr for 
an instantaneous burst assumption. The 3 a upper limit of 8 
Myr corresponds to a minimum progenitor mass of 25 Mq, 
in line with the predictions o f the collapsar mo del for the 
origin of long duration GRBs l|Heger et al.ll2003l ). 

The main limitation is related to spatial resolution (for 
z = 0.1685 one HST/ACS/WFC pixel corresponds to 140 
pc) and arises because we cannot be sure that the GRB 



parent population can be described by a single stellar pop- 
ulation. Thus, nearby clusters of different age may contam- 
inate the photometry, thereby mimicking an extended star 
formation history. Allowing for a more extended star forma- 
tion history, most models still predict a minimum progenitor 
mass of around 25 Mq, but for a small number of models the 
minimum allowed (within 3<r) progenitor mass is ^ 12Mq. 

We discuss ways to obtain tighter age constraints and 
identify the equivalenth width of Ha and near-IR photome- 
try as useful additional diagnostics. 

If exercised on a sufficient number of long duration 
GRBs, the method in this paper could provide information 
on the progenitor mass and metallicity distributions. This 
would provide an independent comparison to estimates ob- 
tained through modeling of the associated supernova. 
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